CO 

~--J  AMRL-TH-65-40 

to 

rH 


O 


RESPONSE  OF  A  LINEAR  DAMPED  DYNAMIC 
SYSTEM  TO  SELECTED  ACCELERATION  INPUTS 


STANLEY  RARRETT 
PETER  R.  PAYNE 

FROST  ENGINEERING  DEVELOPMENT  CORPORATION 


TSpnZCTqF 
HARD  COPY 
MICROFICHE 


rr 

j 


$.  'J.r< 
$.  t.  -  f 


A PHIL  1965 


DDC 


DOC  iRA  £ 


BIOPHYSICS  LABORATORY 
AEROSP  \CF  MEDICAL  RESEARCH  LABORATORIES 
AEROSPACE  MEDICAL  DIVISION 
AIR  FORCE  SYSTEMS  COMMAND 


NOTICES 


When  US  Government  drawings,  specifications,  or  other  data  are  used  for  any 
purpose  other  than  a  definitely  related  Government  procurement  operation,  the 
Government  thereby  incurs  no  responsibility  nor  any  obligation  whatsoever, 
and  the  fact  that  the  Government  may  have  formulated,  famished,  or  in  any 
way  supplied  the  said  drawings,  specifications,  or  other  data,  is  not  to  be 
regarded  by  implication  or  otherwise,  as  in  any  manner  licensing  the  holder 
or  any  other  person  or  corporation,  or  conveying  any  rights  or  permissioi  to 
manufacture,  use,  or  sell  any  patented  invention  that  may  in  any  way  be 
related  thereto. 


Requests  for  copies  of  this  report  should  be  directed  to  either  of  the  addressees 
listed  below,  as  applicable: 

Federal  Government  agencies  and  their  contractors  registered 
with  Defense  Documentation  Center  (DDC): 

DDC 

Cameron  Station 
Alexandria,  Virginia  22314 


Non-DDC  users  (stock  quantities  are  available  for  sale  from): 

Chief,  Input  Section 

Clearinghouse  for  Federal  Scientific  &  Technical  Information  (CFSTI) 

Sills  Building 

5285  Port  Royal  Road 

Springfield,  Virginia  22151 


Change  of  Address 

Organizations  and  individuals  receiving  reports  via  the  Aerospace  Medical 
Research  Laboratories  automatic  mailing  lists  should  submit  the  aadressograph 
plate  stamp  on  the  report  envelope  or  refer  to  the  code  number  when  corre¬ 
sponding  about  change  of  address  or  cancellation. 

Do  not  return  this  copy.  Retain  or  destroy. 


800  -  May  1965  -  448-41-904 


RESPONSE  OF  A  LINEAR  DAMPED  DYNAMIC 
SYSTEM  TO  SELECTED  ACCELERATION  INPUTS 


STANLEY  BARRETT 
PETER  R.  PAYNE 


FOREWORD 


This  study  was  initiated  by  the  Biophysics  Laboratory  of  the  Aero¬ 
space  Medical  Research  Laboratories,  Aerospace  Medical  Division,  Wright- 
Patterson  Air  Force  Base,  Ohio,  with  the  support  of  the  Biodynamics  Section, 
Environmental  Physiology  Branch,  Life  Systems  Division,  National  Aeronauti¬ 
cal  and  Space  Administration,  Manned  Spacecraft  Center,  Houston,  Texas.  The 
research  was  conducted  by  Frost  Engineering  Development  Corporation,  3910 
South  Kalamath  Street,  Englewood,  Colorado,  under  Contract  No.  AF33(657)- 
9514.  Mr.  Peter  R.  Payne  was  the  principal  investigator.  Mr.  James  Brinkley 
of  the  Vibration  and  Impact  Branch  was  the  Contract  Monitor  for  the  Aero¬ 
space  Medical  Research  Laboratories,  while  Mr.  Harris  F.  Scherer  was  the 
NASA  liaison  representative. 

The  information  presented  in  this  report  is  a  summary  and  extension 
of  the  analytical  solutions  for  the  response  of  a  dynamic  system  used  to 
represent  the  mechanical  response  characteristics  of  the  human  body.  The 
solutions  described  within  this  report  were  required  to  supplement  a  research 
program  concerned  with  the  investigation  of  human  body  support  and  restraint 
system  dynamics. 

The  research  presented  in  this  report  was  initiated  in  April  1963  and 
completed  in  June  1963. 

This  report  is  catalogued  by  Frost  Engineering  Development  Corporation 
as  Technical  Report  194-10  and  5s  one  of  a  series  of  reports  generated  in  the 
area  of  human  restraint  and  support  system  dynamics. 

This  technical  report  has  been  reviewed  and  is  approved. 

J.  W.  HEIM,  PhD 
Technical  Director 
Biophysics  Laboratory 


ABSTRACT 


The  general  theory  is  developed  for  the  response  of  a  single  degree  of 
freedom  dynamic  system  to  an  arbitrary  acceleration  forcing  function.  Closed- 
form  solutions  are  obtained  for  a  variety  of  discrete  pulse  shapes  using  the 
method  of  Laplace  transforms  and  the  form  of  the  solutions  indicated  for  oscil¬ 
latory  Inputs  and  semi- Infinite  ramps ,  in  terms  of  complex  Fourier  series. 

A  comparison  of  base  and  mass  excitation  of  the  system  is  included.  In  pre¬ 
viously  published  work  on  this  subject ,  analytical  solutions  are  in  general 
only  given  for  undamped  systems;  an  exception  is  the  response  to  a  sine-wave, 
which  appears  in  many  standard  texts .  The  dynamic  analysis  of  the  human 
body  usually  considers  models  involving  damping,  so  that  in  this  area  there  is 
a  definite  need  for  the  extensions  given. 
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SECTION  1 


INTRODUCTION 


This  report  is  concerned  with  the  behavior  of  a  spring/mass/damper 
system  of  the  type  shown  in  Figure  1.  This  model  has  been  used  in  the  past 
to  simulate  the  dynamic  response  of  parts  of  the  human  body,  in  efforts  to 
assess  human  tolerance  to  acceleration. 

The  response  of  this  system  to  a  continuous  sinusoidal  input  is 
frequently  treated  in  engineering  texts;  however,  the  solution  for  other 
types  of  input  is  usually  confined  to  an  undamped  system.  Mindlin,  et 
al  (Ref.  1)  have  presented  results  showing  the  response  to  a  half-cycle  sine 
pulse,  'lthough  an  analytical  solution  was  given  for  the  zero-damping  case, 
the  curves  for  non-zero  damping  were  plotted  from  analog  computer  results. 

In  Ref.  2  it  is  pointed  out  that  in  the  absence  of  damping  the 
base  excitation  can  be  expressed  in  terms  of  displacement  or  velocity  as 
alternatives  to  acceleration.  Although  this  is  also  possible  when  damping 
exists,  the  resulting  equations  contain  more  than  one  input  derivative  and 
are  generally  too  complicated  to  be  considered  for  engineering  applications. 

He  shall  therefore  confine  our  attention  to  the  response  of  the 
system  when  the  base  of  the  spring  is  accelerated,  but  it  Is  shown  that  the 
results  are  applicable,  by  a  simple  substitution,  to  the  situation  in  which 
the  mass  is  excited  while  the  spring/damper  element  is  grounded. 


SECTION  2 


THE  GENERAL  THEORY  OF  A  LINEAR  SYSTEM 


2 . 1  Comparison  of  Base  and  Mass  Excitations: 

Consider  the  basic  svstem  (Fig,  1)  with  the  excitation  = 
applied  at  the  base.  Then,  if 

•  • 

U  =  acceleration  of  the  mass,' 

W  1*\ 


*8 

=  force  in  the  spring, 

2K  8 

=  force  on  the  damper, 

we  have 

•  1 
mu 

=  fe8 

or  3„ 

=  oar* &  v  2c  8 

since 

09-  -  )_!?_ 

J  rr. 

and 

C  *  -K- 
rr\ 

Let  X 

=  length  of  spring  when 

and  S 

=  deflection  of  spring 

'(1) 


MASS 

(m) 

DAMPER 

(2K) 


Figure  1.  Single  Degree  of 
Freedom  System. 


so  that  &  =  "A 

\ 

iS- 

3 

l 

cC 

n 

• 

8  - 

•  • 

•  •  %  # 

-  3m  + 

•  • 

v 

-  ^  +  9c 

(3) 

Substituting  eq. 

(3)  into  eq.  (2) 

and  rearranging  gives 

*  • 

3,  • 

i  +  2c  8  -f 

0*8 

(4) 

Co..bider  now  the  case  in  which  the  mass  experiences  a  force  excita¬ 
tion  ■P  =  7(t)  ,  with  the  base  of  the  svstem  fixed. 


Corresponding  to  eq.  (2)  we  have,  in  general. 


+  <*i>  +  *2cS  <5> 

Since  the  base  is  fixed, 

Ue  =  constant, 

“<*  at  =  9c  =  ° 

eq.  (3)  becomes 


•  •  •• 


and  substitution  into  eq.  (5)  leads  to 


-IL  *  S  +  2cS  +  (6) 
m 

This  indicates  that  solutions  obtained  for  the  base  excitation  con¬ 
figuration  can  be  applied  to  mass  excitation  problems  simply  by  substituting 
(■Vm  )  for  . 


2.2  General  Solution: 

For  a  linear  system,  we  have  seen  that 


*  *  • 

&  +  2c  &  +  &  =  yc(fc)  (7) 

Using  the  transform  =u*’t  ,  eq.  (7)  becomes 

&"  +  2c  s'  +-  £  *  ('t)  <8> 


where  the  prime  is  used  to  indicate  differentiation  with  respect  to  the  non- 
dimensional  time  parameter T  ,  and  C  = 

This  equation  is  most  conveniently  solved  for  specific  cases  by  the 
use  of  the  operational  calculus.  The  Laplacian  transform  is  defined  by 

QO 

=  f  e  P  x(ar)  dT  o) 

Jo 

using  the  symbol  X(X  )  to  denote  the  transform  of  the  function  x(T  ). 

The  existence  of  the  transform  must  be  verified  when  the  function  to  be  trans¬ 
formed  is  not  common  enough  to  be  tabulated.  Three  sufficiency  conditions  are 
available;  these,  if  satisfied,  guarantee  the  existence  of  (  X  ): 
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(a)  x (*C  )  is  at  least  piecewise  continuous  in  any  interval 

«Tl4*t4,r1  ,  where  *1,  >0 

(b)  <Z'  "Kt)I  is  bounded  near  =  0  for  some  n<  1. 
(de^kr)!  is  bounded  for  large  values  of  *C  ,  for  some 


number  M . 


is  bounded  for  large  values  of  *C  ,  for  some 


In  Laplacian  notation,  eq.  (8)  becomes 


<*)[  t>a  +  2*P  +']  =  +&o('P+'2c)+(€,/)0  (10) 

where  60  =  initial  deflection  at  T  -  0 

and  =  initial  velocity  at  T  =  0. 

We  shall  concern  ourselves  only  with  subcritical  damping,  where 
C  <  1.0.  Thus  we  can  employ  the  substitution 


and  write 


4-  Q.c^>  I 


l-c2,  (11) 

'p1’  4-  *2  C  p  4“  "5  4*  ( l  —  3  ) 

(*p4-  zy  4-  rf  (12) 


so  that  eq.  (13)  yields 

'se(s)  =  ^(ac)  ,  »), 

{p+c)1-  +  V  (• P*2)1  +  '01  (v+zf  +Y)1 

Taking  the  inverse  transformation  and  making  use  of  the  linearity 
of  the  operator, 


b(x)  =  W>~T  c")  1  t  Sn^'T 

C'P*  E)1  +  VX  (p  +  l)1 

■*-($')„•?£''  _ ! _  <u> 

[fp+3)l+V  _ 

In  the  second  term  of  eq.  (13), 

(p4-  C)*1  4-  T)*1  (*p  4-  TT)*  (14) 
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-  giving,  after  inversion. 


C 

V 


S\r\  T 


) 


where 


(f>  =  orcsin  ^ 


For  the  initial  velocity  term, 

(*%*'[**'+  ^  ] 


un0 

V 


-CT 

e  sin  nj«r 


(15) 

(16) 


(17) 


Thus  the  solution  to  eq.  (13)  may  be  written  as 

&  .=  >£(  He") 

When  the  motion  due  to  an  initial  condition  only  is  required,  this 
can  be  obtained  immediately  from  eq.  (18).  Moreover,  for  any  forcing  func¬ 
tion,  *he  initial  condition  transients  are  directly  additive  to  the  response 
due  to  the  forcing  function.  Thus  we  need  only  consider  the  variation  in  eq. 
(18)  due  to  the  particular  forcing  function  in  what  follows. 

We  shall  see  in  Section  3  that  any  acceleration  input  pulse  or  vi¬ 
bration  has  a  characteristic  peak  value  V*  =  \  ( 

Ur'  \ 

It  is  convenient  to  express  &  in  terms  of  Vc*  ,  so  that  eq.  (18)  becomes 


maw 


>-i 


afar) 


l 


+  initial  conditions  solution. 


Initial  conditions  solution: 


+  (&')o  "‘('h* 


(19) 


Consider  the  steady-state  pax t  of  the  solution  given  in  eq.  (19): 


*_  _  *-«  r  *0£)  k-vc-'i  _  ■  wr  l 

yt*  ~  J  X*  l('P*E)l+'9l  J 

The  Laplace  transform  of  eq.  (20)  is  simply 


^(s) 


*(si) 

G>«*)x+,ol 


(20) 


(21) 
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so  that 


I 


(22) 


(ifM-zy+yf 


•A(i)  may  be  regarded  as  the  transform  of  the  "output" 
tem,  with  *^e")  the  transform  of  the  "input";  thus,  subject  to 
initial  conditions,  we  may  define 


of  the  sys- 
quiescent 


w  "  -  ^i) = 

n 

-  as  the  system  transfer  function  for  the  input  yc  ,  with  the  in¬ 
put  applied  to  the  base  of  the  system  as  shown  in  tig.  1.  It  should  be  noted 
that  definition  (23)  is  only  unique  when  the  point  of  application  of  the  in¬ 
put  is  specified. 


The  transfer  function  of  the  system  is  the  Laplace  transform  of  its 
rbeonomic  normal  response  (i.e.,  its  response  as  an  explicit  function  of 
time,  for  quiescent  initial  conditions)  to  a  first  order  unit  impulse. 

From  eq.  (23),  it  follows  that 


^(t)  =  Fi(?)  <24> 

M*)  -  t 

=  •’ [G»(?)  Fj  (-p)3 


Now  the  inverse  transform  of  the  product  of  two  transforms  can  be 
obtained  by  means  of  the  convolution  theorem: 

U-'  [a(^f(-v)]  =  J\(r'-s)jl(s)ds  ) 

=  {  9(s)^(<i-s)ds  j 

where  _ 

9(T)1 

=  [ 

r,(v) 

-  ^(yc’O 

and  S  is  an  arbitrary  "dummy"  variable  of  integration. 
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Thus,  when  an  acceleration  input  function  yc  is  given,  the  re¬ 
sponse  of  a  linear  system  to  y^  may  be  found  by  integration  of  one  of  equa¬ 
tions  (26),  making  use  of  the  known  response  of  the  system  to  a  unit  impulse. 
Although  the  technique  is  based  on  the  Laplace  transform  approach,  it  is  not 
necessary  to  use  the  transform  of  the  input  acceleration  function.  Thi  meth¬ 
od  is  probably  more  useful  than  the  direct  inverse  transformation  pr<  idure 
for  cases  where  the  transformation  is  not  a  rational  function  of  ,  j 

that  the  ratio  c'O/tp  4-  ty-  +■ t)  cannot  be  expanded  into  simple  partial 

fractions.  In  many  cases,  however,  the  integrand  is  so  complicated  that  the 
integration  becomes  very  laborious  to  perform  analytically,  and  numerical 
methods  are  preferable. 

The  approach  just  discussed,  involving  the  convolution  integral, 
may  be  developed  even  further  for  the  special  case  of  a  steady-state  oscil¬ 
latory  input.  The  input  may  then  be  expressed  as  a  complex  Fourier  series: 

°*  Ln-rt/E 

a'  =  «>  •  I c" e  <27) 

n«-  oo 

Where  Cn  is  the  nth  complex  coefficient  in  the  series,  and  Jt,  is  the  frequency 
of  oscillation  of  the  input.  Cn  is  given  by 


C„ 


(28) 


The  complex  exponential  form  of  the  Fourier  series  thus  needs  only 
one  formula  for  all  its  coefficients,  and  is  more  compact  and  easier  to  ma¬ 
nipulate  than  the  more  common  trigonometrical  form.  It  also  provides  a  con¬ 
venient  transition  to  Fourier  integrals  and  Fourier  transforms. 

Re-writing  eq.  (27)  for  and  substituting  into  the  second  of 

equations  (26)  leads  to 


i.'nn.'C 

since  is  independent  of  6 

For  convenience,  let 


f°°  ~  Ln-a(T-s) 

*  J  9(S)  A,  C"e  (29) 

o  n«- oo 

22  in-n/t  f°*  /  \ 

-  Icne  J  9 (6)e  <A*  (30) 


-In  jt.  s 
e 


ds 


(31) 


where  is  thus  the  Fourier  transform  of  the  impulsive  response  function 
sc*)  .  Since  $(*)  is  zero  f or  T <  O  ,  the  limits  of  integration  can  be  0,0® 
or  +  oo  without  affecting  the  value  of  the  integral. 
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Hence  we  can  write  (30)  more  succinctly  a3 


(32) 


For  an  input  which  is  real,  the  output  must  also  be  real;  under 
these  conditions  eq.  (32)  becomes 


0O 

=  C„  A(o)  +  2X  Cn  A(n)cos(nna  +  pn)| 


r\m\ 


where  is  a  phase  angle,  defined  by 


tan 


▼/A 

.  j  9/a  S»n  r>nJl  dT 


I  i,(T)c«njvT:  a*r 


(33) 


(34) 


Since  any  periodic  function  can  be  expressed  as  a  Fourier  series  to 
a  high  degree  of  accuracy,  the  response  of  a  single  degree  of  freedom  system 
to  any  periodic  function  can  be  calculated  from  eq.  (33).  Fourier  series 
representations  of  some  periodic  functions  are  given  in  Table  2  of  Section  3. 

It  should  be  noted  that  the  Fourier  series  expansion  of  a  periodic 
function  which  has  finite  discontinuities  leads  to  finite  (though  not  usually 
significant)  amplitude  overshoots  in  the  region  of  the  discontinuities.  For 
example,  for  a  rectangular  wave,  the  Fourier  series  expansion  gives  amplitude 
overshoots  at  the  discontinuities  which  approach  a  limit  of  about  +  9%  of  the 
amplitude  of  the  wave  as  the  number  of  terms  in  the  expansion  tends  to  in¬ 
finity  (see  Fig.  2).  At  points  on  the  rectangular  wave  between  the  disconti¬ 
nuities  the  difference  is  much  less,  of  course,  decreasing  as  the  number  of 
terms  taken  in  the  expansion  is  increased.  This  behavior  is  due  to  the  fact 
that  the  Fourier  series  expansion  fails  to  converge  uniformly  at  the  discon¬ 
tinuities;  this  is  known  as  the  Gibbs  phenomenon.  It  is  discussed  in  detail 
in  Ref.  3. 


Input  functions  which  are  nonperiodic  (e.g.,  acceleration  pulses) 
may  also  be  treated  by  a  modification  of  the  Fourier  series  expansion,  in  the 
form  of  Fourier  integrals  and  Fourier  transforms.  Generally  speaking,  how¬ 
ever,  the  method  of  Laplace  transforms  usually  proves  to  be  more  convenient 
to  use  for  such  functions. 

Further  discussion  of  the  Fourier  methods  is  beyond  the  scope  of 
this  report.  A  more  detailed  treatment  will  be  found  in,  for  example. 

Ref.  3. 
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Figure  2.  Fourier  Series 


ro|H 


t 

j 


ions  for  Rectangular  Wave. 


SECTION  3 


DISCUSSION  OF  ACCELERATION  INPUTS 


and 


Most  idealized  acceleration  inputs  fall  into  one  of  three  classes: 

(i)  discrete  pulses 

(ii)  oscillatory  wave-forms 

(iii)  semi-infinite  ramps. 


In  this  paper,  the  first  class  has  been  investigated  in  greatest 
detail;  using  the  Laplace  transform  technique,  the  response  of  the  basic  sys¬ 
tem  to  several  different  types  of  pulse  has  been  analyzed.  For  each  pulse, 
a  family  of  curves  is  presented,  showing  the  variation  of  amplification  fac¬ 
tor  with  pulse  duration,  for  a  range  of  values  of  the  damping  coefficient 
ratio.  These  are  plotted  in  Figures  4  through  8.  Amplification  factor  is 
defined  here  as  the  ratio  of  maximum  dynamic  deflection  to  static  deflection. 
The  pulses  considered,  together  with  their  Laplace  transforms,  are  summarized 
in  Table  1  of  this  section;  the  detailed  solutions  will  be  found  in  the 
Appendices . 


Oscillatory  inputs  may  also  be  solved  using  Laplace  transforms; 
however,  a  more  convenient  approach  in  general  is  by  way  of  the  Fourier  se¬ 
ries  expansion  of  the  input  function.  This  method  is  discussed  in  detail  in 
Section  2  of  this  paper.  The  Fourier  series  representations  of  some  periodic 
functions  are  presented  in  Table  2. 

Semi-infinite  ramps  represent  acceleration  inputs  which  start  from 
an  initial  value  of  zero,  and  rise  in  some  fashion  to  a  constant  acceleration 
level.  Rise-time  is  the  most  important  parameter  for  this  type  of  input. 
Table  3  shows  some  of  the  ramps  commonly  encountered,  together  with  their 
Lap lace  transforms . 
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TABLE  1:  DISCRETE  PULSES. 


TYPE 


DIRAC  IMPULSE 


RECTANGULAR  PULSE 


WALE-SINE  PULSE 


VERSED -SINE  PULSE 


TRIANGULAR  PULSE 


FUNCTION 

L  -  TRANSFORM 

yc  =  (t*i  &) 

l 

*(•0“  Av 

yc  =  V*  (o<'t<a) 
=  O  (T  >a) 

^9;')=Yt"(i-e-ap) 
T>  V  ' 

9"  =  Yc"  sin  At 

(o<x<  2) 

Jl' 

=  0  (T>£) 

*feC) 

=  YfV  (l+e-^) 

/ 

*(»:)  « 
y;  (^r(i-e_ap) 

(0<T<a) 

=  0  (T>a.) 

2  HpM  w 

9c"»^*T  (o<T:<a) 

=^'(2.a- T)  (a<T<lj 

a.  ' 

=  0  (^>Zo) 

*(9,*)  = 

TABLE  2  :  REPRESENTATION  OF  SOME  PERIODIC  FUNCTIONS 
.  BY  FOURIER  SERIES. 
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b  +T 


TABLE  3:  SEMI-INFINITE  RAMPS. 


TYPE 


EXPONENTIAL  FRONT 
•  • 

«c 


FUNCTION 

L- TRANSFORM 

y„  =  Yc  (t>o) 

^(bt)=l 

yc^  1  (o<b<t^ 

®VC  (t>tR) 

^('O  = 

1-  (l-e-pt) 
t„-px  V  > 

U  s  4  (  1  -  COS  TT  fc  ) 

C  ^ 

(o<  fc<  tR) 

*(V  = 

Y.  RL3'(l+e^) 

=  YC  (t>t„) 

ae=Yc(l-e-at) 

(t  >o) 
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APPENDIX  I 


DIRAC  IMPULSE  FUNCTION 


From  Section  3,  the  equation  for  the  Dirac  impulse  function  is 

simply 

ij*  =  &V  (35) 

with  the  Laplace  transform 


■e£(4c'')  = 


In  (19),  (replacing  by  *Ar) 


(36) 


_ l _ 

_  -eTT 

•  Qy  6  _  _C _  Skicv  OjT 

by  Av  **) 


(37) 


(38) 


Equation  (38)  is  plotted  in  Figure  3  for  several  val^^s  of  c  .  It 
is  interesting  to  note  that  in  the  initial  phases  of  the  motion  it  is  ade¬ 
quately  described  by  a  system  consisting  of  a  mass  and  damper: 


i.e. 


,  -2cT 

I-  e 


) 


(39) 


As  the  deflection  increases,  the  spring  force  becomes  important 
enough  to  cause  a  divergence  from  the  simple  damper  result,  and  returns  the 
system  to  equilibrium,  after  passing  through  a  maximum  deflection. 


We  determine  the  raaximun  deflection  by  differentiating  eq.  (38) 
with  respect  to  *t  ,  then  equating  to  zero. 


-CT  _  -C*C 

c  cos  nr)T  —  JL  e  Sm^T  «=  o 


/.  I  -  £_  7?<r*=  o 

Solving  fort*,  the  time  at  which  maximum  deflection  occurs. 


(40) 

gives 


=  *  arctan  *0 

» 

at*s'»n  7) 

7)  -c 

(Note  that  ^"L^- 

- ►  1  as  C  - 1 

«c«- 

^  IT  -  _ 

^  as  C  ' 

►  O) 

t> 
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Substituting  for*t*in  eq.  (38), 

t*"  femax 

US'  Av 


Mote  that 


and 


US  Av 


I  as 


^  are  plotted  in  Pig.  4. 

US'  AV 


US 

US' 


When  the  initial  deflection  is  not  zero,  eq.  (38)  may  be  wri 

-Z'Z 


us 


!l_.  =  e'  ^  yt +[c  t  +  4.]  sin  yr  j 

ing  the  substitution  *f  =  ^  for  brevity. 


US  AV 

where  &0  -  initial  spring  deflection. 

This  can  be  simplified  to 

us1^ 


(A  AV  ^ 

where  ©  =s  Q1XSH1 


-cT  f - - 

= _ +  2c  $  +■  I  sin('»)'r  + 


JTTziFPi-’ 

Differentiating  with  respect  to  *T  , 


d 

J  'S7  2  c  1  + 1  5one  CTcxss//ncr  +  ©) 

dT 

|  us  Av  J 

r  w 

-cT 


=  e 


—  c  e  s\n(yi  +  e)  j 

J'S1+'2el  +  i  sln^T  +  e  +  y) 


where  sin  (fl  =  HQ  as  before. 


(41) 


tten 

(4:) 


(43) 


(44) 
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Equating  to  zero,  this  simplifies  to 


1)T  4-  ©  4-  Of) 

Sin  (y)  <T*  4-  e) 


Also 


*  mr 

=  .sin  (  hit  —  <fi) 

=  ±fY) 

HTT  -  8  -  9* 


(45) 


(46) 

(47) 


Substituting  eq.  (46)  and  (47)  into  eq.  (43)  gives 


WLlsaai  _  e  ^(V"e)jtl+2cS  + 


where 


(a»  AV 

sin  0 
sin  <jp 


Jl?  (o«0«l))(48) 


i'S7  +2zt  +1  .  _ - 

*9 


The  response  of  an  undamped  system  to  an  impulsive 
for  various  values  of  the  initial  parameter  ^  ,  is  given  in 


velocity  change, 
Figure  5. 
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AMPLIFICATION  FACTOR,  ^Av** 


Figure  3.  Response  of  a  Damped  Single  Degree  or 
Freedom  System  to  an  Impulsive  Velocity  Change. 
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DAMPING  COEFFICIENT  RATIO  c 


Figure  4.  Variation  of  Amplification  Factor  and  lime  of  Maximum 
Deflection  with  Damping  Coefficient  Ratio,  for  an  Impulsively 
excited  Linear  dystem. 
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Figure  b.  Influence  of  Initial  Deflection  Rarameter 
on  Response  of  the  System  to  an  Impulse. 
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APPENDIX  II 


RECTANGULAR  PULSE 


Prom  Section  3, 


In  eq.  (19) , 


-{ 


>c 

=  v; 

ii  O 

(o^n:  a.) 

(49) 

^  ( a") 

=  VC#  / 
?  V 

'1  —  e”°p) 

(50) 

s 

=  *r' . 

l-c"- 

(51) 

Vc" 


-C'l 


-  |  4-  e~CL  sin/noT  —  arcJton  V\ 

*)  -2  ' 


\  -f 


-c('T-a)  % 


*) 


sm 


^/Y)pT-a]  -  arcban  H^T-a) 


and 


where  H(T-q)  denotes  the  Heaviside  unit  function: 

H  (T-cx)  =  O  ,  fT<  a. 


J 


T  >  a 


S  =  i 

a 


-cT 
_  Sm 


jn  ^cr)T  —  arctun  ^0  ^ 

-c 


^or^a  (52) 


&  _  J 5  e  sin/ljT  —  arcban  *9  ) 

yc"  V  \  v  -c  / 

^sin^pT  —  a. J  -  crcck-an  5. )  JorT>a 


d_ 

dT 


-e 

Differentiating, 

/JL\  =  * 

k  yc*J  Vc* 


(53) 


-CT  ,  / 

— jL  e  Sinf^T  -  anctran 
+  e~cT  ocs(^'t  —  arctan 

=  e  C  sin  yr 


5) 

-c  • 


(54) 
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and  for  the  second  time-regime, 


-cT  -c(*-a) 

—  sitl  ‘pT  —  ? _  sin  TJ  (t-o^ 


>  a!)  (55) 


We  can  find  T  ,  the  time  at  which  occurs  by  setting  eq.  (55) 

equal  to  zero,  ana  solving  for  time.  The  value  of  Sm*M  is  then  found  by  sub- 
stitu  ing  *T  =  *t*  into  the  equations  for  S. 


( i;  for  T<  at  ,  setting 


6>yc" 


=  0  gives 


Sin  =  0 

Of)1!*  =  nit 


In  (53),  _  I  4-  e  77  sin  ( mr  -  arcbon  2.  ) 

TT  V  ~z 


*  -  ^ 


and  for  n  =  1,  we  have 


—  C  IT 


^max  =  l  +■  €  5 


(T  <  a) 


V  " 

Tt  s 

(Note  that  for  zero  damping,  we  get  the  familiar  result  "*2^"  =  2,  indicating 
100%  overshoot  of  the  system.) 

(ii)  For  *r  >  O.  t  eq.  (55)  leads  to 

e”cT  |  sin  or)T*  -  e?a  S»n  ^(T*-  cx)|  —  O 


and,  ignoring  the  trivial  result  of  *t  =  oo  ,  we  get 

*  OCX  (  *  A  ,  *\ 

sin  r)  T  =  e  P  SmnpT  Cos  TJO.  —  CoS'TjT  sin  TpO,C 


Dividing  through  by  Sin 

i  =  e'4 


i#  * 

SmnpT  CcrS  'tJO.  —  CcS'TjT  sin  T)a 

* 

»  # 

:d  jeers  Of)  a  —  cot  'pT  sin'TJOLj 


-COL 


cxytT)Tw  = 


cos  nr)  a.  —  © 

s»v»  t?0. 


whence 


^  =  1  arc  tan/  sin  T)a  ) 

'*)  (  cos^QOl  -  c-cd  j 


In  (53), 


& 


max 


arefctxn 


(59) 


The  response  of  a  damped  single  degree  of  freedom  system  to  a  rec¬ 
tangular-pulse  acceleration  input  is  plotted  in  Figure  6,  as 
versus  the  non-dimensional  pulse  duration,  c 


Using  Figure  4,  it  is  possible  to  calculate  the  initial  slope  of 
the  Dlot  of  jwwi  vs.  (Ora.  for  any  arbitrary  pulse,  providing  a  useful  check 

on  the  solutionj  this  check  was  applied  to  the  case  of  the  rectangular  pulse. 
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NON-DIMENSIONAL  PULSE  DURATION,  uf  a 

Figure  6.  Response  of  tne  System  to  a  Rectangular  Pulse  (Zero  Initial  Conditions). 


APPENDIX  III 


HALF-SINE  PULSE 


The  equation  of  the  pulse,  from  Section  3,  is 


^  =  Vc  sin  (o  £  ~) 

,  O  (*>7i)  (60> 


with  the  Laplace  transform 


-PIT 


=vt 


*) 

'p1,  4-  -A.a 


Substituting  eq.  (61)  into  eq.  (19)  gives 


yc“ 


«C 


-i)-n.(l  +  e  ~"~  ) 

(p  Xil^pS-'lcp  +  l) 


In  terms  of  tne  Heaviside  unit  function, 


vc" 


=  =  I—  ( at  -  V',) 

J(i -,n>)'*,+4c\n>  |  -a 


-CT 

x  f= _  S»n 

V 


(61) 


(62  > 


-*(*-£) 


The  phase  angles  are  O.CJL 


'V*.  =  (Vrctran 


%  =  avcbati  ^Z1  (rT«^TLn)(6l4) 


We  must  consider  separate  solutions  for  the  two  time  regimes 


a)  <r<  v, 


Because  of  the  Heaviside  unit  function,  eq.  (63)  becomes 


—  -  =  4-e 

Vc"  J(»-A')a+  4c'a'  l-a  ") 


(65) 


Note  that  for  7L  -  0,  this  simplifies  to 


I  —  SI? 


sm  aT  -  a  sm 


ioTj 


since  “i*t— »»0  and  ►IT  as  C  — — ►  0.  This  agrees  with  the  zero-damping  solu¬ 
tion  for  established  in,  for  example.  Ref.  4.  Differentiating  (66), 

v/  y1-^+4cvn  v  v  ] 


tfhere 


'V'j  =  arxXxxn (^P)  —  atesm  If)  ^ 


(67) 


—  ^  * 

*1  ,  the  time  at  which  6  achieves  its  maximum  value,  is  the  solu¬ 
tion  of  the  equation  obtained  by  setting  (67)  to  zero: 


-ft- _ ^  casing*- ’fr)  +  e _  O  (68) 

J(i-a')V4sV™(  ^  ^  V  J] 

An  explicit  solution  of  eq.  (68)  in  *C  ,  if  it  exists,  is  certainly 
not  simple  to  find,  except  for  the  case  of  2  =0,  when  eq.  (68)  becomes 

* 

OOPS  Aft  -ockT  =  0 

so  that 

ccri  SI  T  -  CJCTS  T 

.AT*  -  0.T m  ±  T* 


whence,  for  T3  =  1, 


T  - 


0.TT 

a  ±  I 


Since^C  is  necessarily  positive,  and  A  may  be  less  than  unity,  we  must  take 
the  positive  sign  in  the  denominator  of  eq.  (69). 

For  non-zero  values  of  c  ,  is  found  most  directly  by  a  numeri¬ 
cal-graphical  approach. 

(ii)T  >  %. 

For  this  case,  since  (T-Vn.  )  >  0,  the  value  of  the  Heaviside 
unit  function  is  unity,  and  eq.  (68)  becomes 


on 


Noting  that 


=  F====  \ 
J(l-Aa)l  +  42'jta  ( 


+ 


_L_  sm(.a*c  -t,) 
-a 


sin^T-t) 


♦^r!KA[''SH) 


(70) 


we  may  simplify  eq.  (70)  to 

&  -ft- _ e~c*J 

<"  1  +  c5Va  (71) 

Differentiating  eq.  (71)  with  respect  to*T  gives 

S'  -ft  e~c<  f  -C  sin /tit:  -  ^ 

yc*  J(l-^V4cW  V  \ 

-ce^sin('f  +1)  <*s('’)'1  '^) 

+■  Sin  (<!) ' T  -  X  -  T)  VjO  ]  (72  ) 

Now  CBS  X  —  c  sin  X  =  s^X+arctan^) 

=  sin(x  +  V,) 

by  the  definition  of  ^  in  eq.  (67). 


.*.  (72)  becomes  z_  . 

I' _ -n.  e  Jsinfr'C  -%+%) 

yc‘  “  oo-^y+^e’jv1,  ?>  1 

sinA)T— ,yi+')*a_71^j»)^ 

V  (73) 
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Equating  to  zero,  we  have  the  equation  for  : 


sin(nr)T# -/V«X  4-*^)  -+*  e  — *0, 


(74) 


ignoring  the  trivial  solution  of  €  =0. 

We  can  again  make  the  useful  check  of  inspecting  the  zero  damping 
case:  for  6=0,  eq.  (71)  bec'mes 


s 

V 


=  -n.  £  sin  (t:  -  it)  +  s»n  -  tv  —  } 

1 -  ( 

=  _  |  _  sin  <T  -  s»n(eT-7^n)| 


(75) 


The  solution  for  this  case  given  in  Ref.  4  may  easily  be  rearranged  to  give 
eq.  (75), 

Setting  E  =  0  in  eq.  (75)  gives 

Sin  *  —  TV  +  4-  sin  (t*  —  TT  -t-  -  O 

COS  <1*  4-  COS  (  T  *  -  ^Sl)  ~  ° 


or 


T*  =  rvn  ±  (T*  -"^dn.) 


'2't*  =  hit  ±  ^ 


_n_ 


Finally,  for  n  =  1, 


<T*  = 


TT 

Z 


_l_  TT_ 

Z  -0. 


(76) 


Substituting  for*T  into  eq.  (72)  leads  to 

^mox  _  ~  il  i  COS  IT 

x  1  2X 


v.1 


>--o>  l  2ji 
-2  SI  ar&  IT 


I  -.TL1  <2jX 

£ 

which  agrees  with  the  solution  for  given  in  Ref.  4. 


(7/) 


We  need  to  solve  eq.  (74)  explicitly  for  T  . 


For  brevity,  write 


a  = 


so  that  eq.  ( 7*0  becomes 

sin^Of)T*-A)  +e  SinMf)T*  -3) 


Z*/Jl  . 


(78) 


=  o 


sin  nr)fC  cos  A  —  ccs  ijT  sin  A 


j^n  HQ  X  *  Cess'S  —  CCS  T)T*  sin  6  j  =  O 


Sin  T)T*|  CCS  A  +  €  ^  COS 


=  CCS  nf)T*|  sin  A “ft‘sin^Bj 


.'.  tan  nr)T 


CTJ/ 

sin  A  +  e  sin  3 


whence 


ecs  A  4-  e  ^  ccs3 
ss  I  one  ban  f  sin  A+e  sin  3  } 

7)  1  ccs3  4-  e^’^cos'B  J 


or,  substituting  for  A  and  B, 


T  «s  1  onefcon  sin(l^x  (79) 

7>  «»(V**)-*-  e'^cw^vv?**’) 

Maximum  deflection  may  now  be  derived  by  substituting  for  *1  =  T*in  eq.  (71). 
The  response  of  the  system  to  a  half-cycle  sine  pulse  is  plotted  in  Fig.  7, 
for  a  range  of  values  of  damping  coefficient. 
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NON-DIMENSIONAL  PULSE  DURATION,  or  a 

Figure  7.  Response  of  the  System  to  a  Half-sine  Pulse  (Zero  Initial  Conditions). 


APPENDIX  IV 


VERSED-SINE  PULSE 


From  Section  3, 


ll 


£fl  - 
2  L 


era  2txT 
a  . 


(oc'C^o.) 

(  t>a) 


(80) 


^£(yt")  _  t  !(■ 


-ap\ 

i-e)  (si) 


Using  eq.  (19),  and  writing  as  D,  for  convenience. 


J> _ LaT’f 

I  -  e-“p 

ril-eT0*)  ) 

Vc"  ~  z  | 

/F>(VX-*»'2cp  1  j 

fp*+-J>'][pV'iep*i] ) 

where 


and 


A. 

vc" 


A 

vcrt 


say 


Si 

V' 


I  +  £■ 


V 


•  *r‘i 

l-e-ap  } 

Z  |  l)  ) 

(  -CT  . 

hT-^l  +  [] 

•  }  e  Sin 

2  |  g) 

l'  '}  \ 

^  s'n  fo[T-c] 

1 

-e-^)  ) 

^  1  (fx ^xXi>a -*-0  i 

=  -L  _  |  cxrs(j)<Z-’fS) 


(82) 


(83) 


(84) 
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The  phase  angles  are 


arcfccm 


>  =  osrdcoirx 


•&) 

(sW) 


=  -T&  —  arck-an 


\  D*1,  +  2.7:  -  l  J 


The  complete  general  solution  of  Eq.  (82)  is  thus 


e -  Sin(/Y)ct -Vf)  +  J_ 

2t)  V  j  J  z 


fi 


-c  -0.1 


~  lL  - -Ljccs/yi-^+el!  sln^T-*,) 

-  gots|i>[t  -a]  h(t  -  a) 


e  S<n 


We  now  separate  the  general  equation  into  the  twc  time  regimes  ,*T<  a.  and  *t>a 


(b<T<  a 


h(<t -a)  =  o 


so  that  eq.  (86)  reduces  to 

r  -cT 


i_  -  e 

X,"  1y)  v  ' 


-  ^  ^  cos/dt  -*■  elx  sln/^T ■-•)$ 

z  1  v  v  v  yj 
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Differentiating, 


_  5-Diin(l)T 

J^TWl  v  ' 

+  e'5Tf-|.  J 


53  e _  sm  -rjT 

*T> 

+  ^  ^ 

J(lM)a+-4z'3’-1  ( 

k  O  when  T  =T  and  ^  **  >rnax 


A  graphical  method  of  solution  was  used  to  solve  f  or  *1  *  ,  and  Smftx  found  by 
substituting^  =T  •  into  eq.  (87). 


00*1  >  a 

For  this  case,  h(t-»)  =  1,  and  eq.  (86)  becomes 

_S_  =  £lZ  s'.nO)<c-V;)-ee*-T  ^W'9[T-a]->;) 

yc"  ii)  v 

_  ^X  f e~CT  sin^T-^ 

+v^ tt  HtH"  4)} 

(89) 


Differentiating  V  with  respect  to  *1  and  simplifying  the  result  gives  the 
following  result : 
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&max 


is  found  by  substituting  equations 
eq.  (89). 


(90) 


and  (91)  into 


The  response  of  the  system  to  a  versed-sine  pulse  is  plotted  in  Figure  8. 
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Figure  8.  Kesponse  of  the  System  to  a  Versed-sin.  Pulse  (Zero  Initial  Conditions). 


